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5h \ Abstract. We review a recent construction of an explicit analytic series representation for 

symmetric polynomials which up to a groundstatc factor are eigenfunctions of Calogero- 
Sutherland type models. We also indicate a generalisation of this result to polynomials 
which give the eigenfunctions of so-called 'deformed' Calogero-Sutherland type models. 
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The purpose of this paper is to give a pedagogical review of a recent construction of an explicit 
analytic series representation for symmetric polynomials which up to a groundstate factor are 
eigenfunctions of Calogero-Sutherland type models. In special cases, this construction has been 
previously studied in [HJ [12] , while a detailed account of the more general results presented 



in this paper is in preparation [8]. 



To establish a context for our construction we begin this introduction by briefly discussing 



a 

, quantum many-body models of Calogero-Sutherland type in general and highlighting some of 

the distinguishing features of those models which have eigenfunctions given by polynomials. By 
reviewing Sutherland's original solution method for the Calogero model [23] we proceed to recall 
that these polynomials have a triangular structure and to discuss its importance when explicitly 
constructing them. We then sketch the main steps in our solution method and compare it to 
Sutherland's. An outline for the remainder of the paper finally concludes the introduction. 



1.1 Quantum many-body models of Calogero-Sutherland type 

A quantum many-body model of Calogero-Sutherland type is for some potential functions V 
and W defined by the Schrodinger operator 




where N refers to the number of particles present in the system, and Xj to their positions. As first 
observed by Calogero [3] and Sutherland [23] in two special cases there exist certain choices of 

*This paper is a contribution to the Proceedings of the Workshop on Geometric Aspects of Integ- 
rable Systems (July 17-19, 2006, University of Coimbra, Portugal). The full collection is available at 
|http: / / www.emis.de/journals / SIGMA/Coimbra2006.html | 
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these potential functions for which the resulting Schrodinger operator can be exactly solved. For 
many of these choices this is due to the fact that its eigenfunctions are given by natural many- 
variable generalisations of polynomials solving a second order ordinary differential equation. 
This includes for example the original models of Calogero and Sutherland, whose eigenfunctions 
respectively correspond to the Hermite polynomials and the ordinary monomials p n = x n . We 
also mention Olshanetsky and Perelomov's [19] root system generalisations of these models, 
where the Legendre, Gegenbauer and Jacobi polynomials similarly appear; see [T| I26j. 

All of these models share a number of remarkable properties: their square integrable eigen- 
functions, labelled by partitions A = (Ai, . . . , Ajv), i.e. integers Aj such that Ai > • • ■ > Aat > 0, 
are all of the form 

*aOi, • • -,xn) = ^o(zi, ■ ■ ■ ,x N )P x {z(xi), . . .,z(x N )) (2) 

with particular symmetric polynomials P\ and a ground state which always is of the product 
form 

TV 

V (xi, ...,x N ) = Y[ Mxj) Y[{z{xk) - z{xj)) K , (3) 

3=1 j<k 

where the function z is fixed by the choice of potential function V, and tpo is the ground state 
of the one-body model obtained by setting iV = 1 in (pQ). The corresponding eigenvalues are 
in addition of a very simple form and can be written down explicitly. In the remainder of 
this paper we will refer to the polynomials P\ as reduced eigenfunctions of the corresponding 
Schrodinger operator ([1]) and our aim is to explain an explicit analytic series representation for 
them, obtained in [8]. 

1.2 Triangular structures and Sutherland's solution method 

If we conjugate the Schrodinger operator ([1]) by its groundstate (J3]), and subtract the corre- 
sponding eigenvalue Eq, we obtain the differential operator 

j=l J 3=1 J J 

which has the symmetric polynomials P\ as eigenfunctions. It was observed already by Suther- 
land [23] that a key property in their construction is that this differential operator can be 
consistently restricted to certain finite dimensional subspaces of the symmetric polynomials, 
on which it can be represented by a finite dimensional triangular matrix. This reduces the 
problem of constructing the reduced eigenfunctions of the Schrodinger operator (pQ) to that of 
diagonalising a finite dimensional triangular matrix. 

To make this more precise we now present a slight modification of Sutherland's original 
argument for the so-called Calogero model, defined by the Schrodinger operator ([1]) for V{x) = 
x 2 and W(x,y) = 2k(k — l)(x — y)~ 2 with k > 0, i.e., 

N F> 2 N 1 

j=l ° X j j=l j<k [X 3 Xk > 

Note that we without loss of generality have set the harmonic oscillator frequency w > to 1: 
it can be introduced by scaling Xj — > ^fooxj and H — > luH. It was shown by Calogero [3] that 
this Schrodinger operator has eigenfunctions of the form ([2]), with groundstate 

TV 

vPo(xi, . . . ,x N ) = n e" (1/2 H l[(x k - x 3 r (6) 

j=l j<k 
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corresponding to the eigenvalue 
Eq = N(1 + k(N-1)), 

and where the symmetric polynomials P\ are natural many-variable generalisations of the Her- 
mite polynomials. This implies that these symmetric polynomials are eigenfunctions of the 
differential operator 

H N :=m Q \H N -E )m 

dx ) + 2 X ° dx i ^ ^ x i ~ Xk ( dx J dxk ) ' 

The idea is now to compute the action of this differential operator on the so-called monomial 
symmetric polynomials m\, for each partition A = (Ai, . . . ,Xn) defined by 



m\(xi, x N ) = ^ x i PW 



X N ■ 



where the sum extends over all distinct permutations P of the parts Xj of the partition A. In 
the discussion below we will on occasion refer to monomials m n parametrised by integer vectors 
n G Nq which are not partitions. Such a monomial is then defined by the equality m n = m p ( n ), 
where p(n) refers to the unique partition obtained by permuting the parts rij of n. Note that 
as A runs through all partitions of length at most N the monomials m\ form a linear basis for 
the space of symmetric polynomials in N variables. Using the fact that 

^- + 2x-^] x n = 2nx n - n(n - l)x n ~ 2 , 



dx dx 
as well as the identity 



\dx~~d^) {xnym+ynxm) = {n ~ m) S x^y^-mix^y^+y^x™- 1 ), 



x — y\dx dy J 

valid for all x, y € R and n,m G No such that n > m, it is straightforward to verify that 

N L(A 3 -A fc )/2J 

H N m x = 2\X\m x -^2*j(Xj - l)m X -2e 3 ^ ~ X k)m X -^ + i)e J+ {u-i)e k 

j=l j<k v=\ 

+ 2k^2 hmx- ej -e k , (8) 
j<k 

where |A| =A + -- - + Aat, [n/2\ denotes the integer part of n/2, and ej are the natural basis 
elements in 7j N defined by (ej)j. = Sj^. It is important to note at this point that the right hand 
side of this expression involves terms which in general are not parametrised by partitions, e.g., if 
A = (3, 2, 2) then A — 2e\ = (1, 2, 2) which is not a partition. However, since mpn) = m x, for all 
permutations P of N objects, we can remedy this problem by collecting all terms corresponding 
to the same monomial m\. Once this is done we find that the action of the differential opera- 
tor ((7J) on the monomials m\ is triangular, in the sense that if two partitions \i = (/xi, . . . ,/j,n) 
and A = (Ai, . . . , Ajv) are ordered according to the partial ordering 



fj, < X 4^ fii-\ h Hj < Ai H h Xj, V j = 1, . . . , N, 



4 



M. Hallnas 



then Hm\ is a linear combination of m\ and monomials with \x < A and |/i| < |A| — 2, i.e., 



for some coefficients c\^, and where the sum is over partitions [i < A such that < |A| — 2. 
This means that when constructing the reduced eigenfunctions of the Calogero model we can 
restrict the differential operator ([7]) to a subspace of the symmetric polynomials spanned by 
monomials m^, where fj, < A for some fixed partition A. On this subspace the differential 
operator (|7j) can indeed be represented by a finite dimensional triangular matrix, with off- 
diagonal elements ca m , and where its diagonal elements 2|A| give the eigenvalues for the reduced 
eigenfunctions of the Calogero model, which correspond to the eigenvectors of this matrix. 
There remains then to actually compute the matrix elements c\^, i.e., to collect all terms in (JSj) 
corresponding to the same monomial m^. It seems however that this problem does not have a 
simple solution, which in turn implies that the reduced eigenfunctions of the Calogero model 
do not have a simple series representation in terms of monomial symmetric polynomials. The 
situation is similar for the other models discussed above (see e.g. [8]), and as far as we know 
also for other simple bases of the space of symmetric polynomials, such as elementary, complete 
homogeneous and power sum symmetric polynomials; see e.g. [T7] for their definition. 

1.3 A sketch of our solution method 

To obtain our explicit analytic series representation for the reduced eigenfunctions of Calogero- 
Sutherland type models with polynomial eigenfunctions we use a construction which differs from 
the one discussed above in two important aspects: first, we express them in terms of a particular 
set of symmetric polynomials /„, n G Ta N , on which the action of the differential operator (j4]) 
is simpler than on the symmetric monomials; second, we avoid the problem of computing the 
matrix elements analogous to the in ([9]) by using an overcomplete set of these polynomials, 
parametrised not only by partitions but by a larger set of integer vectors in Z Ar . One could of 
course apply this latter change to the symmetric monomials and the discussion in the previous 
section. Note, however, that the expression ([8]) for the action of the differential operator (J7J) 
on the symmetric monomials is valid only for partitions, and that a formula valid for arbitrary 
integer vectors in Nq would be more involved. 

To simplify notation we will here, and in the remainder of the paper, let x = (x\, . . . ,xn) and 
y = (yi, ■ ■ ■ ,Vn) be two sets of independent variables. For an arbitrary integer vector n £ T, N 
we will furthermore use the notation X — X • • • X j^- ■ and similarly for y. We now define the 
set of symmetric polynomials f n , n G Z , through the expansion of their generating function 



valid for \y^\ > ••• > \y±\ > m&Xj(\xj\). Although the expansion unavoidably generates terms 
parametrised by integer vectors which are not partitions, we prove in Section [23] that a basis for 
the space of symmetric polynomials is formed by those f n which are parametrised by partitions 
alone. The reason that we use precisely these symmetric polynomials is that for each Schrodinger 
operator ([!]) there exists an identity 




(9) 




(10) 




(11) 
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where CV is a constant, H N is obtained from by a simple shift in its parameters (see [8]), 
and the function F is given by 



where i/jq is the groundstate of the one-body model obtained by setting N = 1 in Hj^ . Note 
that if the groundstate factors are removed and the variables Zj = z(xj) and Wj = z(yk) are 
introduced we essentially recover the generating function for the symmetric polynomials f n . 
This relation will later enable us to obtain the action of the differential operator on the f n 
in a straightforward manner. As is then shown, this action is simple enough to be inverted 
explicitly, thus yielding our explicit analytic series representation for the reduced eigenfunctions 
of Calogero-Sutherland type models with polynomial eigenfunctions. 

In the literature there exist various other approaches to the construction of these reduced 
eigenfunctions. In a recent paper Lassalle and Schlosser [16] obtained two explicit analytic 
series expansions for the Jack polynomials, the reduced eigenfunctions of the Sutherland model, 
by inverting their so called Pieri formula. For very particular partitions or a low number of 
variables explicit analytic expansions have also been obtained by other methods; see e.g. [T7] . 
In addition, various representations of a combinatorial nature are known for the Jack, as well 
as certain other related many- variable polynomials [51 [91 \T7\ [27] . We also mention the recent 
separation-of-variables approach to the Sutherland model due to Kuznetsov, Mangazeev and 
Sklyanin [10] , which also relies on the identity (jlip . This list of previous results reflects only 
those which we have found to be most closely related to ours. For a more comprehensive 
discussion we refer to [8]. 

1.4 An outline for the remainder of the paper 

We continue in Section 2 to give a more detailed account of our solution method by applying it 
to the particular case of the Calogero model. In Section 3 we then discuss generalisations of this 
result to other Calogero-Sutherland type models with polynomial eigenfunctions and also to 
the 'deformed' Calogero-Sutherland type models recently introduced and studied by Chalykh, 
Feigin, Sergeev and Veselov; see [U [211 E2] and references therein. 

2 A first example: eigenfunctions of the Calogero model 

In this section we provide a detailed account of our solution method by applying it to the 
Calogero model, defined by the Schrodinger operator ([5]). Apart from the proof of completeness 
these results were all obtained in [7J. 

We begin by formulating our main result: an explicit analytic series representation for the 
reduced eigenfunctions of the Calogero model in terms of the symmetric polynomials f n . In 
doing so we make use of a few notational conventions which we now introduce. In contrast 
to the introduction we will here use the following partial ordering ordering of integer vectors 
m, n G Z N : 

m -<n ■£> rrij + • • • + m/v < rij + ■ ■ ■ + n^, Vj = 1, . . . , N. 
To simplify certain formulae we associate to each n £ Z N the shifted integer vector 




K 



11 



+ 



(n 



+ 
i ' 




nj = rij +k(N+ 1 -j). 
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For each integer vector n G % N we define the Kronecker delta 

N 



i=i 

We also recall the notation ej for the standard basis in "L N , i.e., (ej) k = 5j k . We are now ready 
to state the main result of this section. 

Theorem 1. For an arbitrary integer vector n G 7L N let 

Pn = fn + '^2 U «( m )/m, (12) 
m 

where the sum is over integer vectors m G Z N such that 

m <n and \m\ < \n\ — 2, 
and the coefficients 

oo.. oo / s \ s / r \ 

«»M=Ei^ E-E E M"»+E^ 11^ --"-E^ . 

s=l ' ii<fci j s <fc a ^i,...,^ s =0 V r=l / r=l V £=1 / 

where we use the shorthand notation 

gjk(v, m) = 2k(k - l)u(l - 5 jk ) - m+(m+ + l)5 u0 5 jk (13) 

and 

E v jk = (1 - u)ej + (1 + u)e k . 

Then P n is a reduced eigenfunction of the Schrodinger operator ([5]) corresponding to the eigen- 
value 

E n = 2\n\+E , E = N{1 + k(N-1)). 

Moreover, as A runs through all partitions of length at most N the P\ form a basis for the space 
of symmetric polynomials in N variables. 

Remark 1. It is important to note that the series defining the coefficients u n (m) terminate 
after a finite number of terms, and thus are well-defined. This is a direct consequence of the 
definition of the Kronecker-delta 5 n (m) and the fact that the equations 



n — m = E 



l'r 

r=l 



only have a finite number of solutions v = (yi,..., u s ) for fixed n,m G 



7i\ 



To prove the theorem we proceed in three steps: we begin by deriving the identity (jlip for 
the Schrodinger operator ([5|); we then prove the first part of the theorem, that the functions P n 
are reduced eigenfunctions of the Schrodinger operator ([5]); and finally, we prove that a basis for 
its eigenspace is given by those eigenfunctions which are parametrised by a partition of length 
at most N. 
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2.1 The identity and a model with different masses 



Rather than proving the identity (jlip for the Schrodinger operator (|5|) by a direct computation 
we obtain it here as a consequence of a more general result which has the interpretation of 
providing the exact groundstate of a generalisation of the Calogero model where the particles 
are allowed to have different masses. We will, however, not stress this interpretation but rather 
use the result to derive various other identities, of which (jlip is the one of main interest for the 
discussion which follows. 

Proposition 1. For a given set of real non-zero parameters m = (mi, . . . ,mj^) and variables 
X = (Xi,...,Xat) let 



M i d 2 N i 

i'=l 3 3 j=l j<k 3 



(14) 



and let 



j=l j<k 

We then have that 

W$o = £o$o (15) 

with the constant 

/ JV \ 2 Af 

\j=i J j=i 

Moreover, if all rrij are positive and <£(> is square integrable then fi defines a self-adjoint operator 
bounded from below by £q and with groundstate &q. 

Proof. We prove the statement by establishing that the differential operator (|14p is factorisable 
according to 

3=1 3 

with 

Qf = ±d Xj + Vj, V 3 = ^dx 3 $o- 

Note that is the formal adjoint of Qj . The identity (|15|) then follows from the fact that 
Qj&o — for all j. If all rrij are positive then this factorisation shows that 7i defines a unique 
self-adjoint operator via the Friedrichs extension which is bounded from below by £q (see e.g. 
Theorem X.23 in [20J) and with <I>o as ground state. 
Observing that 

Vj(Xi, . . .,Xjv) = -rrijXj + KS^m ] m k 
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it is straightforward to deduce that 

3=1 3 

with remainder term 

Upon symmetrising the double sum and using the identity 

1 1 1 

+ T-Z ~ TTTT 7TT + 



(Xk — Xj)(Xj - Xi) (Xi — Xk)(Xk - Xj) (Xj — Xi)(Xi — Xk) 
it is readily verified that 

N \ 2 M 



( N \ z N 
k^j \j=l J 3=1 



So. 



We note that by setting all rrij = 1 we obtain as a direct consequence of the proposition 
that © indeed is the groundstate of the Calogero model. On the other hand, setting M = 2N, 
rrij = 1 and rriN+j = — 1 for j = 1, . . . , N we see that Ti splits into a difference of two Schrodinger 
operators ([5]) and that we obtain the corresponding identity (|1 1 j) with hL ^ = H^. 

Corollary 1. With 

n Uivk- VjY 

3 1 j,k 

we have that 

H N (x)F(x,y) = (H N (y) + C N )F(x,y), (16) 
where the constant 

C N = 2(1 - k)N. 

It is interesting to observe that Proposition [1] implies a number of additional identities. We 
can for example choose to take different number of variables Xj and yk- This leads to an identity 
involving two Schrodinger operators Hn and Hm with different number of variables N and M. 
We may also set some of the parameters rrij to either 1/k or —1/k while still preserving the 
property that Ti splits into a difference of two differential operators, which in this case will 
define so-called 'deformed' Calogero-Sutherland type models; see Section [3721 These additional 
identities are further discussed in [8]. 

2.2 Construction of reduced eigenfunctions 

We proceed to prove the first part of the statement in Theorem [TJ that the symmetric polyno- 
mials P n , defined by (fT2j) . are reduced eigenfunctions of the Schrodinger operator ([5]). We begin 
by computing the action of the differential operator ([7]) on the symmetric polynomials f n . 
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Lemma 1. For each n 6 7L N we have that 

N oo 
H N fn = E n fn ~ ~ " 2 )fn~2e 3 + 2k{k - 1) £ £ uf n _ {1 _ v)e ._ {1+1/)ek (17) 

j=l j<k v=l 

with 

E n = E n - E = 2\n\. (18) 

Proof. We first note that the function F in Corollary [1] and the generating function for the 
symmetric polynomials f n are related as follows: 

f» x n (i- ^ 

F(x, y) = *„(*) II ^-i(%) y 7 K(7V+1 " j) J<k ) Vk 

\j=i ) n(i-- 

with VI/o the groundstate © of the Schrodinger operator ([5]). The identity (|16p in Corollary [1] 
together with definitions ([7]) and (I10|) . of respectively the differential operator i/^r and the 
symmetric polynomials f n , therefore imply that 

(H N f n (xj) y~ n+ = Y, fn(x) {H N + C N - E ) y~ n+ , (19) 

where 

N 



h n = -jj H N fl ik,-i(Vj) 



n ^o,-i{yj) j=1 

3=1 

N d 2 



We now expand the interaction term in a geometric series 



(20) 



v- 



E 



(»-») 2 frTyr 1 



which in the region \yn\ > \yN—i\ > ••• > \yi\ is valid for all j < k. It is now straightforward 
to compute the right hand side of (fT9|) . and by comparing coefficients of y~ n+ on both sides of 
the resulting equation we obtain (fT7|) with 



N 

E n = Y( 2n t - 1) + Cjv - E . 

3=1 

As a simple computation shows, this indeed coincides with (118p . and the statement is thereby 
proved. ■ 

Remark 2. At this point it is interesting to compare the action of the differential operator ((7J) 
on the monomials m\, given by ([8]), and on the polynomials f n , as just obtained. We note, 
in particular, that the simpler structure of the latter arise from the fact that it is essentially 
equivalent to the action of the differential operator Hn, defined by (|20|) . on the powers y~ n+ 
and the fact that the 'interaction' terms of this operator does not contain any derivatives, in 
contrast to the differential operator (J7|). 
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It is clear from (|17|) that the action of the differential operator (|7|) on the symmetric poly- 
nomials f n has a triangular structure, in the sense that Hf n is a linear combination of f n and 
symmetric polynomials f m with m -< n and \m\ < \n\ — 2. This suggests that to each n £ 
corresponds a reduced eigenfunction P n of the form (|12|) with eigenvalue £" n . Inserting this 
ansatz into (fT7|) and introducing u n (n) = 1 we obtain 

/ iv 



-Hat-P„ = £■„/„ + ^ E m u n (m) - + l)m+ii n (m + 2ej 



m \ j=i 

oo \ 

vu n {ra + (1 - u)ej + (1 + v)e k ) 

fmt 

j<k v=\ / 

where the sum is over integer vectors m 6 1i N such that m -< n and |m| < |n| — 2. We therefore 
conclude that the validity of the Schrddinger equation Hj\r^ n = E n ty n follows from the recursion 
relation 

oo 

2(\n\ - \m\)u n (m) = ^ XI 9jk(v, m)u n {m + Ej k ), 

with the coefficients gj k (y\m) defined by f)13j) and where we used the fact that E n — E m = 
2{\n\ — \m\). We now proceed to solve this recursion relation. Suppressing the argument m we 
rewrite it in the form 

Un = <5n + Ru n , 

where the operator R is defined by 

1 



(Ru n )(m) 



2(|ra| — |m|) 



^2 ^29jk(v, m)u n (m + E u jk ). 



j< k u=0 



Observe that this expression is well-defined since \n\ — \m\ ^ for all applicable m and the sum 
truncates after a finite number of terms: u n (m) is by definition non-zero only if m ^ n. The 
solution of this latter equation is therefore 



u n = (l-Ry 1 5 n = J2 RS5 n 



s=0 



where the expansion into a geometric series is well-defined since it only contains a finite number 
of non-zero terms, as will become apparent below. Using the definition of the operator R as 
well as the defining properties of the Kronecker delta 5 n we deduce that 



(R s 5 n )(m) = £ £ 



j B <k s ia,=0 



2(|ra| — |m|) 



EE 

ji<ki ux=o 2( In 



E E 

js-l<k s -l Vs-1=Q 
S 



gjs-iks-^s-i-im + E"*) 



2[\n\ - \m + E: 



jsk s 



9j 1 k 1 [vi;rn+ E 2k t 



jiki 



5 » m+ E E £ 



E -E E *. »+E E ^ n 



r=l 

9j r k r [vr;n- £ E j ek 



31<kl js<ks U lt ...,U 3 =0 



r=l 



=1 2 n 
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By finally observing that 

r 



4r 



we obtain our explicit analytic series representation (|12p for the reduced eigenfunctions of the 
Calogero model. 

2.3 Completeness of the reduced eigenfunctions 

There remains only to prove that the reduced eigenfunctions just obtained provide a basis for 
the space of symmetric polynomials, i.e., that they span the eigenspace of the differential opera- 
tor ([7|). We obtain this last part of Theorem [1] by exploiting the relation between the symmetric 
polynomials f n and the so-called 'modified complete' symmetric polynomials g\, defined through 
the expansion of their generating function 

r ^5 A (x)m A (y" 1 ), 



n(i- a ) » 

j,k\ VkJ 

valid for min^. \y^\ > maxj(|xj|), and where the summation extends over all partitions of length 
at most TV. It is well known that the g\ are homogeneous symmetric polynomials of degree |A|, 
and also that as A runs through all partitions of length at most N they form a basis for the 
space of symmetric polynomials in N variables; see e.g. Section VI. 10 in [T7]. We mention that 
these first properties can be directly inferred from their generating function, whereas the fact 
that they span the space of symmetric polynomials is a consequence of the equivalence between 
the expansion by which they are defined and the fact that they are dual to the to the monomial 
symmetric polynomials m\ in a particular inner product; see e.g. Statement 10.4 in |17j . 
By comparing the generating functions for the f n and the g\ we find that 



j<k 

Assuming \y^\ > ■ ■ ■ > \y±\ and expanding each term in the product in a power series we rewrite 
the right hand side as follows: 

K \ ~n+ E Pjk(ej-e k ) 
j<k 



nGN^ j<kp jk =0 

where we have taken p(n) to denote the unique partition obtained by reordering the parts nj 
of n. This means that 

OO , v 

/« = II E (-^'ij^E^ler^))' (21) 

j<k Pjk =0 ^ kJ 

Many of the properties of the 'modified complete' symmetric polynomials g\ for this reason 
carry over to the f n . In particular, observing that 



n + ^pjkiej - e k ) 

j<k 
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for all integers pjk we conclude that they are homogeneous symmetric polynomials of degree n. 
We also see that they are non-zero only if n y 0. Now suppose A is a partition. It is then clear 
that 

m := A + ^Pjk(ej - e k ) < \ 

j<k 

and furthermore that also p(m) ^ A: we obtain the partition p{m) from m by some permutation 
of its parts rrij, and since by definition fx\ > ■ ■ ■ > ^ijy for any partition fi = (/zi, . . . , ^n) we 
have that p(m) < m < \. We therefore conclude that 

h = g\ + J2 M ^9» ( 22 ) 

for some coefficients M\^, and where the sum is over partitions fi -< A. As indicated in this 
expression we let M = (Mx^) denote the transition matrix, defined by the equality f\ = 
M\^g^, from the fx to the g^. Given a partition A it follows from (|22p that it can be 

consistently restricted to the partitions \i such that ^ A. With rows and columns ordered in 
descending order this restricted transition matrix is upper triangular with l's on the diagonal. 
Hence, it can be inverted. Since the inverse of an upper triangular matrix is upper triangular 
we obtain that 

5a = /a + ^(M" 1 ) A/1 / / „ 

where the sum is over partitions [i -< A. We have thereby proved the following: 

Proposition 2. The functions f n are non-zero only if n y 0. In that case, f n is a homogeneous 
symmetric polynomial of degree \n\. Moreover, as A runs through all partitions of length at 
most N the f\ form a basis for the space of symmetric polynomials in N variables. 

The same line of reasoning can now be applied to the reduced eigenfunctions 
-Px = fx + ^2u x (m)f m 

which are parametrised by partitions A. Recall that the sum is over integer vectors m -< A. 
Using formula (|2ip and following the subsequent discussion we find that 

for some coefficients b\^, and where the sum now extends only over partitions fi -< A. Applying 
the arguments leading up to Proposition [2] we conclude that this expression can be inverted 
to yield each g\ as a linear combination of the with /u H A. Hence, as A runs through all 
partitions of length at most N the P\ form a basis for the space of symmetric polynomials in N 
variables. This concludes the proof of Theorem [TJ 

3 Generalisations to other models 

In this section we indicate how the results on the Calogero model obtained in the previous section 
can be generalised to similar models with polynomial eigenfunctions, including not only models 
of Calogero-Sutherland type but also the 'deformed' Calogero-Sutherland models introduced 
and studied by Chalykh et.al.; see [H \21\ [22] and references therein. A detailed account of these 
results is in preparation [8]. 
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3.1 Calogero— Sutherland models with polynomial eigenfunctions 

When the number of particles are set to one in the Calogero model it reduces to the very well 
known harmonic oscillator, which features eigenfunctions given by the classical Hermite poly- 
nomials. This is only one special case of the following well known and more general statement: 
to each complete sequence of polynomials {p n : n G No}, obeying a second order ordinary 
differential equation, there is a corresponding Schrodinger operator 

h = ~& + V{x) (23) 
with a particular potential function V such that its eigenfunctions are of the form 

tp n (x) = ipo(x)p n (z(x)) 

for some functions tpo and z. This can be verified by first observing that such a set of polynomials 
are eigenfunctions of a differential operator 

h = a(Z) d? + P(Z) d-Z> 

where 

ct(z) = a2Z 2 + ol\z + ao and (3(z) = f3\z + /?o 

for some coefficients aj and (3j. Now introducing the variable x = x(z) as a solution of the 
differential equation 

x'(z) 



V a ( z ) 

and defining the function by 



w(z{x)) _../ _ a>' -2(3 



i> (x) = e~ w(z[x>) , w 



4a 



it is straightforward to verify that a Schrodinger operator h of the form (|23p is obtained by 
conjugation of h by the function ip$ and changing the independent variable to x, or to be more 
precise, 

- d 2 
h = -V'O^V'o 1 = ~q^2 + v ( x )> 

with potential function 

t r ( \ , , w (20 - a')(20 - 3a') 1 „ 1 a , 
V(x) = v(z(x)), v = ^ ^- '- + -a" - -0. 

To illustrate this general discussion we have listed the particular values of the coefficients aj 
and /3j, as well as associated functions ipo and z, which correspond to the classical orthogonal 
polynomials (of Hermite, Laguerre and Jacobi) and the generalised Bessel polynomials in Table[TJ 
In fact, we can by simple translations and rescalings always reduce to one of these four cases. 

We mention that transformations of differential equations of the type described above are fre- 
quently used in the theory of ordinary differential equations of second order; see e.g. Section 1.8 
in Szego's classical book [25J on orthogonal polynomials. It is interesting to note that this 
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Pn(z) 


a(z) 


(3{z) 


ipo(x) 




H n (z) 


1 


-2z 


e 1 




(Hermite) 










r(«) 


z 


(l -\- ] — z 


X^G ^ ^ 2 




(Laguerre) 










Pn (z) 


1 - z l 


b - a - (a + b + 2)z 


sm a+i/z - cos 6+i/ ^ - 


COS z 


(Jacobi) 










y n (z; 1 - 2a, 26) 


z 2 


26+ (1 - 2a)z 


exp(— 6e _x — as) 


e* 


(gen. Bessel) 











Table 1. The particular values of coefficients aj and [3j. as well as associated functions ipo and z, corre- 
sponding to the classical orthogonal polynomials (of Hermite, Laguerre and Jacobi) and the generalised 
Bessel polynomials. 

transformation has a simple and direct generalisation to many variables. Setting the interaction 
potential 



W(x,y) 



a{z(x)) + a{z(y)) 
(z{x) - z(y)) 2 



and keeping the potential function V as given above one verifies that the Schrbdinger opera- 
tor dU) after a conjugation by the function \^o, as defined in ©, and a change of independent 
variables from the Xi to the Zi, as defined above, is transformed into the differential operator 



II 



N 



d 2 N d 1 



N 



i=i 



j<k z i 



a{z 



Following the discussion in Section 1.1 it is straightforward to verify that the action of this 
differential operator on the monomial symmetric polynomials is triangular in the very same 
ordering as in the case of the Calogero model. This means that the Schrodinger operator (pj for 
these choices of potential functions V and W, up to degeneracies in its spectrum, has a complete 
set of reduced eigenfunctions given by symmetric polynomials. We mention that this unifying 
point of view on Calogero-Sutherland type models with polynomial eigenfunctions seems to have 
been little used in the literature, with the notable exception of Gomez-Ullate, Gonzalez-Lopez 
and Rodriguez [6] who, among other things, used this point of view to obtain the spectrum of 
all these models. 

In [8] we show that our construction of an explicit series representation for the reduced 
eigenfunctions, presented in the previous section for the Calogero model, goes through virtually 
unchanged for all these models. In particular, we generalise Theorem Q] to the following: 

Theorem 2. For n £ Z, N , the reduced eigenfunctions of the Schrodinger operator ([I]) are 
formally given by 

Pn = fn+'^2 u n{m)f m , 
m 

where the sum is over integer vectors m G 7L N such that 



m -< n and \m\ < |n| + deg(a) — 2, 
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and the coefficients 

q=r 



U 



»m=EE-E E EM™+E*K II 

i=l J"l<fcl ji<fe;Pl,...,Pi=0^,-,n=l \ *=1 / r=l h [n-Y E Vq l q 

\ < * In fcn 



jqkq 

q=r 



where we use the shorthand notation 

N 



b n (m) = E n -E m , E n = - ^2{a 2 n j (n j - 1) + (ft + 2k(N - j))nj), 

3=1 

gjkip, v\ m) = (1 - Sj k )K(K - l)a p (2v - p) 

- 5jk5 v im+(5 p0 a (m+ + 1) + <5pi(ai(mt + « + 1) - /5b)) 



and 



£^ = (l-i/) ej - + (l-p + i/)e fc . 

if b n (m) /or a// integer vectors m G smc/i i/iai m -< n and |m| < |n| + deg(a) — 2 i/jen 
P n is a we// defined symmetric polynomial. Moreover, if this is the case for all integer vectors 
n £ suc/t iftoi n = A /or some partition A of length at most N then the corresponding P\ 
form a linear basis for the space of symmetric polynomials in N variables. 

Remark 3. It is important to note that the condition b n (m) ^ 0, m -< n and \m\ < \n\ + 
deg(a) — 2, is essential in order for the coefficients u n (m) to be well defined. For generic choices 
of the parameter k and the polynomials a and (3 it is satisfied for all n 6 Z, N ; see [8] for a further 
discussion of this point. 

Remark 4. At this point it is interesting to enquire whether our basis for the reduced eigen- 
functions of the Schrodinger operator (pQ) , as stated in Theorem [21 in applicable cases stand 
in a simple relation to the generalised hypergeometric polynomials of Lassalle |13| [T4"l [T5] and 
MacDonald [18], defined by expansions in Jack polynomials. In the case a = —z 2 and = —z, 
corresponding to the Jack polynomials themselves, one can show that they in fact coincide and 
it seems natural to expect this to be true also for the generalised Jacobi polynomials. Since 
the generalised Hermite and Laguerre polynomials are limiting cases of the generalised Jacobi 
polynomials (see e.g. [l]) this would imply the equivalence also in these two cases. If established, 
this result would in these cases imply a natural orthogonality for the reduced eigenfunctions in 
Theorem [2j At this point these statements are only conjectures and we hope to return to them 
elsewhere. 



3.2 Deformed Calogero— Sutherland models 

There exist an interesting deformation of Calogero-Sutherland type models [4"l I2H [22] . defined 
by the following class of differential operators in two sets of variables x = (x\, . . . ,xn) and 
x = (x±, . . . , Xj^): 



id 2 \ ( <9 2 



K - 1 



+ k(k- l)^2w(xj,x k ) + (1 - K)^W(xj,x K ) H W(xj,x K ), 

j<k j,K J<K 



16 



M. Hallnas 



where the potential function V is obtained from V by a simple parameter shift; see [8]. They 
provide a natural generalisation of the models discussed in the previous section, in that they 
also have polynomial eigenfunctions. To be more precise, they have eigenfunctions which can 
be labelled by partitions A = (Ai, A2, • • •) such that Aat + i < N (see [22]), and are of the form 

*a(Z) S) = Vo(x, x)P\ (z(xi), z(x N ),z(xx), . . . , z(xft)) , 

where the function is given by 



N N 



n (z( Xk ) - z{xj)Y n iz(x K ) - z{xj)fi« 

, = i 7-i U(z{x K ) - z{xj)) 



3=1 J=l 

and the P\ are polynomials in the variables Zj = z(xj) and zj = z(xj). They are however no 
longer symmetric under permutations of all variables but only under permutations restricted to 
the Xj or the xj. In addition, they obey the condition 

d d \ 

— + k— )P X = 

dzj dzj ) 

on the hyperplanes Zj = zj, for all j = 1, . . . , N and J = 1, . . . , N. The corresponding algebra 
of polynomials has been extensively studied by Sergeev and Veselov [2"H |2"2"]. 

In [8] we also construct explicit series representations for the reduced eigenfunctions P\ of 
these 'deformed' Calogero-Sutherland type models. The construction is analogous to the one 
discussed in previous sections, with the difference that the reduced eigenfunctions now are 
expressed in a set of polynomials f n ,h, (n,n) E 1j N+n , defined through the expansion of their 
generating function 

nfi-sVnfi-^V 7 " nfi-f^Wi ^ 



i<fc\ VkJ J<K\ VKJ j,K\ VKJ J } k\ Vk 

El I i-r/-, x i\ K T-r ( + XJ 



ni 1 -^) n(i-^) n(i 

VK ) )k\ Vk) }jc\ VK 

fn,n(x,x)y- n y- h , 

valid for \y^\ > > \vn\ > •■■ > \yi\ > maxj t j(\xj\,\xj\). We mention finally that the 

number of variables Xj and xj as well as yj and yj may be chosen differently in the definition of 
these polynomials f n ,ni thus allowing for a number of series representations to be obtained for 
the same reduced eigenfunction, an aspect of our construction which is further discussed in [8]. 
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